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Nuclear pore complexes (NPCs) are very selective filters that monitor the transport be- 
tween the cytoplasm and the nucleoplasm. Two models have been suggested for the plug of 
the NPC. They are (i) it is a reversible hydrogel or (ii) it is a polymer brush. We propose 
a mesoscopic model for the transport of a protein through the plug, that is general enough 
to cover both. The protein stretches the plug and creates a local deformation. The bubble 
so created (prtoein+deformation) executes random walk in the plug. We find that for faster 
relaxation of the gel, the diffusion of the bubble is greater. Further, on using parameters 
appropriate for the brush, we find that the diffusion coefficient is much lower. Hence the gel 
model seems to be more likely explanation for the workings of the plug. 

The nuclear envelope in all eucaryotes is perforated with nuclear pores . Each pore 

has a selective filter, referred to as the nuclear pore complex (NPC). The NPC is a self-assembled, 
eightfold symmetric ringlike structure consisting of eight copies each of 30-50 different proteins, 
connecting the inner and outer nuclear membranes. It regulates the import and export traffic of 
proteins and has two distinct modes of transport: passive and facilitated. Passive transport is non- 
specific and takes place by ordinary diffusion. Colloidal gold particles with radii up to 4 nm, and 
generic proteins up to 50 kDa in mass, pass efficiently through the NPC in this way [7[. In contrast, 
facilitated translocation allows the passage of objects as large as several megadaltons. Proteins 
having a short amino acid sequence known as nuclear localization signal (NLS) form a complex with 
transportin (a transporter protein rich in hydrophobic units) and are transported in this mode. 
The transport requires specific interactions between the translocating species and constituents of 
the NPC and consequently is highly selective. Gold particles of up to 32-36 nm in diameter are 
able to pass through some NPC if they are coated with nucleoplasmin-importin complexes Q]. 
This suggests that the interaction of the protein-transportin complex with NPC is essential for 
transport. Passage of proteins through the NPC has attracted considerable experimental and 



theoretical attention 
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According to Ribbeck and Gorlich (RG) ^ the central plug 



of the nuclear pore is made of long diblock copolymers rich in hydrophobic phenylalanine-glycine 
(FG) units forming a meshwork. Only the macromolecules which can form hydrophobic contacts 
with the FG units are incorporated into this network and get transported. In an interesting paper, 



Bickel and Bruinsma (BE) [10(] point out that such a model would lead to a lower rate of diffusion. 
BE suggest that the central plug is a reversible polymer gel in a poor solvent and a protein in 
it experiences an extra noise (they call it "chemical noise") arising from the fluctuations of the 
FG contacts and this extra noise enhances the diffusion of the protein within the NPC. Single 
molecule fluorescence microscopy by Yang, Gelles and Musser jj| shows that the protein executes 
random walk inside the central core of the NPC. An alternate model for the plug suggests that it 
is not a gel but a polymer brush 1^. Surprisingly, there have been experimental support for both 
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gel and brush models. In interesting experiments, Prey et al. Il3l . \l4 . llSH have shown that the 
nucleoporins form a hydrogel in vitro, offering support to the model of RG and EE. On the other 
hand a beautiful study by Lim et al. [l^ has found that the proetins, when grafted to a surface 
behave like an un-cross-linked brush. Also, it is known that the interaction between the transportin 
and the FG residues is not just hydrophobic, but involves hydrogen bonding, electrostatic and van 
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18l | and that there are extremely hydrophilic portions in between the 



der Waals interactions 
FG units. 

In the following we study a minimalistic model for the transport in the NPC. Our model is 
quite general, and would be applicable whether the plug is a gel or a brush. The actual values of 
the parameters in the model would depend on whether it is a gel or a brush. We find that the it 
is possible for the protein to diffuse rapidly within a reversible gel, while the diffusion would be 
much slower within the brush. 

We take the pore complex to be infinitely long (end-effects neglected) and shall adopt a con- 
tinuum description for the plug. We use x to denote position along the direction of the axis of 
the NPC. To make the problem one dimensional, we imagine the cross section of the pore to be 
a square, with width Ly in the Y and height Lz{= Ly) in the Z directions. We shall assume 
that the particle has a length 2i?o in the X direction, causes a distortion of height ~ a and fills 
the pore fully in the Y direction. Further, to simplify the analysis, we assume that periodic 
boundary conditions are imposed in this direction, with a period Ly. With these, the problem is 
reduced to two dimensions {X and Z). The size of the distortion needed to create a cavity to 
accommodate the particle will be our important variable. Let ipix) denote the height of the cavity 
in the Z-direction at the position x. The simplest possible expression for the energy of distortion 
Edis would have terms quadratic in 4>{x). Thus E^is = (1/2) j°°^dx {d(i){x) /dx)^ + k(f){xY^ . 
The energy of interaction Em = -Eb + f^^dxQ{x - R) ^{k + AA;)/2{0(x) - a}^ - {k/2)(t){xf^ 
of the particle at R with the plug can compensate for the distortion energy. Q{y) with y = x — R 
is a function that determines the interaction between the particle and the plug. We assume 
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Q{y) to be a symmetric function of y, having maximum value Q(0) = 1. Further, we assume 
Q(ibcxD) = 0. When the particle enters the plug it would have to break hydrophobic contacts that 
may be there, and the energy expenditure for that may be met by formation of new contacts of the 
nucleoporins with the particle. All these together is represented by a single constant Ei,. Thus the 
energy of the system is, to within an additive constant, E[(f), R] = E^is + Eint- We refer to this as 
the particle on a string model for the transport, as the above expression is identical to the energy 
of a stretched string, which is displaced from its equilibrium position by the particle, to form a 
bubble. We assume overdamped, Langevin dynamics for the string and the particle, given by 

5E[<j),R] 



dR 



C 
Cp 

^{x, R,t) is 2x1 column vector with elements (j){x, t) and R. (" is the friction coefficient for the 
string and Cp for the particle. We shall use ^ to denote the transpose. F{x, t) = (/(x, t),g{t))^ where 
f{x, t) and g{t) are Gaussian random forces with mean zero and (/(xi, fi)/(x2, ^2)) = 2(kBT5{xi — 
X2)6{ti — 12) and {g{ti)g{t2)) = 2(pkBT5(ti — ^2)- We define /? = l/ksT and denote the time of 
relaxation of the long wave length oscillations of the string hy t = (/k. A model with similar 
structure has been studied in the context of DNA replication by Bhattacharjee [l9|. 

It is convenient to work with dimensionless quantities, and use the original symbols them- 
selves to denote them. Thus we change to /3E E, f3Eh Eh, x/Rq x, (/)(x)/i?o (Pi^)^ 
a/Ro a, R/Ro R, (3kRl k, (3AkRl Ak, a ^ aPRo, t/pCRo ^ t and (p/CRo Cp- 
The minimum energy configuration for the system obeys the two equations "^^^^^^ = and 
^^1^ = 0. These lead to - k^{x) - (AA;((/.(y) - a) - ka)Q{x - R) = and 

dxQ'ix - R){{k + AA;) {0(x) - a}' - kcj) (x)) = 0. The prime, ' is used to denote differ- 
entiation with respect to x. As the shape of the distortion would obviously depend on where the 
particle is located, it is clear that the distortion would have the form (pciv) where y = x — R. In 
terms of y, the first of the two equations become cr^-^l^ — ktpdv) — {Ak{cj){y) — a) — ka)Q{y) = 0. 
The second equation is easily satisfied by having (f)c{—y) = 4>c{y)- Once 4>c satisfying these 
conditions is found, it may be used to get the minimum energy, Ec = E[(j)c{x — R),R\. The 
minimum energy configuration with the particle stationary at R^, and with the center of the 
distortion coinciding with R^ shall be referred to as having the bubble at R[, and this may be 
specified by the function ^'c(x — Ri,,Rb) = {4>c{x — Rb),Ri))^. Movement of the bubble as a 
whole is described as change of 4'c(x — Rb,Rb) by change of Rb- It is clear that this motion 
is translational as the ener gy of the system is unchanged by this movement. This causes the 
existence of a "zero mode" \2(t\. At a finite temperature, there would be fluctuations which will 
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cause the bubble to execute random motion. It is to be noted that the distortion and the particle 
can fluctuate in opposite directions and therefore Rb and R follow different dynamics. Follow- 
ing the methods of instanton theory [2^, that has been used in a similar context 
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22 , we 



write the state of the system at any time as the configuration with the bubble located at Rb 

plus fluctuations about this conflguration, which we expand in terms of the normal modes about 

00 

the bubble at Rb- Thus, we write = ^'c(x - Rbit),Rb(t)) + X] Ci(t)^i{x - Rbit)), 

1=1 

where — Rb) = {(j)i{x — Rb),Ji)^ are normal modes, the equation for which will be defined 

later. We now substitute the above expansion into the Langevin equation and expand the RHS 

00 . 

up to first order in Ci{t). The result is - Rb{t)),l)^ + E Ci{t)'^i{x - Rb{t)) = 

1=1 

Ci{t) \L'^i{x - Rb{t)) - Rb^'iix - Rb{t))\ + F{x,t)where <f is a 2 x 2 diagonal matrix with 
1=1 ^ 

diagonal elements 1 and Cp. The operator L is defined by 

L^iiy) = 

{-ad^dy^ + k + AkQ{y)} Uv) " Q'iv) {AA;(0,(y) - a) - ak] Ji 

- IZ Q'iy') {^My') - a) - ak} Uy') + Ji IZo dyQ"{y) {{k + Ak)/2 (My) - - imMyf} 

We take "^i to obey the eigenvalue equation L^i = Xi'ifi. It is convenient to introduce an inner- 
product (^'/l^ij = ((/>/|<Ah) + CpJlJh, where = I-oo'^^^ii^)^hi^)- With this defini- 
tion of the inner product, L is a hermitian operator on the space spanned by ^i. Further, 
d/dRb{(t>c{x — Rb),Rb)'^ = {~4>'c{^ ~ Rb)AV is an eigen function of L with an eigenvalue zero, 
as may be easily proved by differentiating ( ^fif'f ^ ) = and ( ^^fe^ I = 0, 
with respect to Rb and writing the results in the matrix form. This is the zero mode of the system. 
After normalizing, this mode may be written as ^'o = (— (/>^(x — Rb), l)V\/ci with c = Cp + {4''c\4>'c)- 
Note that in the spirit of instanton approach the sum over I in the expansion of ^{x,R,t) does 
not include this mode. Further, the fact that Q{y) is a symmetric function, means that one can 
classify the eigenfunctions based upon the symmetry or antisymmetry of (j)i{y). If is 
symmetric, then "^i has the simple form ^'^ = (x^O)^, where Xl is a symmetric eigenfunction of 
the operator (—ct-k^ + k + AkQ{y)]. Further, for this case, Xi = ei, where £/ is the eigenvalue 



9?r 

associated with Xh a fact we will use later. 

On taking inner product of the equation involving Ci{t) with ^o(x — Rb{t)) we get Rb{t) = h{t) + 
{Rb{t)/c) E Ci{t) (^(^'^{(pi^ , where h{t) = {g{t) — (^cp'^ \ f{t)'^)/c. To get the simplest approximation 
for the diffusion coefficient for the bubble, we neglect the "string-phonon-particle scattering" (SPP) 
term (the last term) of the above equation involving Rb{t) and get Rb{t) = h{t). The diffusion 
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coefficient is then found to be Dq = ^ /g dti Jq dt2 {h{ti)\h{t2)) = c~^. One can write Z^o as a 

harmonic mean Dq = DpDg {Dp + Dg)^^ , with. Dp = Cp^, Dg = (^(p'^ \ (j)^ ■ It is interesting 

that within this approximation, Dq is always less than Dp implying that the diffusion coefficient 

within the plug is less than outside, as expected [8]. To get a better approximation, one has to 

include the last term of the equation involving Rb{t). For this, we have to find Ci{t). On taking 

inner-product of the equation involving Ci{t) with ^i{x — Rh^t)), and neglecting SPP scattering, 

we get Ci{t) = —\iCi(t) + {^i I F{t)), which can be solved with the condition Ci{—oo) = to get 

^li^) — /-oo I -^(^i)) ^^P(~-^/(^ ~ should now solve for Rb{t) and then calculate 

the diffusion coefficient. However, this is difficult, as the noise is multiplicative. We therefore adopt 

a simple minded approach in which Rb{t) on rhs is replaced by h(t) to get the following equation 

Rb{t) = h{t) (1 + B{t)/c), where B{t) = YZi too (^^ I ^(^i)) exp(-A;(t - h)) | (^z). If 

one further assumes that h{t) and B{t) are uncorrelated then the diffusion coefficient becomes 

D = Dq + Jo dti Jq dt2 {h{ti)h{t2)) {B{ti)B{t2)). This is easy to calculate since the noises 

are delta function correlated and it becomes, D = (l + -B/c^) , where B = {B{t)B{t)) = 

Yliyto {p^c \ 'PiJ /-^z- As 02 {x) is an even function of x, only 0; that are symmetric will contribute 

in this sum. But as we saw above, they are identical with xu and for them A; = We 

can therefore write B = \xi^ l^i- It is clear that one can expand the sum over I 

-1 



to include all I. Thus we get D = Dg ( 1 + 



G 



0;')/c2). G = [-a^ + k + l^kQ{y) 



is the Green's operator. Note that in the Schrodinger operator —cr-^ + k + AkQ{y), the term 
AkQ{y) vanishes at infinity. It is straightforward to calculate the diffusion coefficient for simple 
models of Q{y)- We define an enhancement factor r = D/Dq = ^1 + (^cp^ G cf)^"^ /c^^ • We 
now perform calculations for a simple model with Q{y) = 9{y — 1) — 9{y + 1). Further, we put 
AA; = so that the calculation becomes easy. Then one can find (pc easily and obtain Ec = 
—Eb + a'^^/kas\l\h.{^^k/a) exp(— y^/c/cj), where Eb is the dimensionless value of the binding energy 
(= Eb(3). If Ec < 0, then the entry into the NPC has no activation energy (however, exit would be 

G 



activated). The Green's function is G{y,yi) = ( y^k/Aa) exp(— -v/^ \y — yi\) so that 



^/k/aj [l - exp{-2y^k/a) - 2y^k/a + 4k/aj /8. With these, we get Dq = 2a/6 with 6 

i exp(2yf ) (^-l+exp(2^)-2yf +4f ) 



2( (-l+cxp(2^))a2 7^-2M^+2Cpexp{2y^) 



{2aCp + a\Vka-exp{-2^J^){2k + Vka))) and r = 1+ 

We now estimate the diffusion coefficient of the protein, assuming that the plug is a gel. As 
experimental information on the properties of the gel is scarce, the numbers used are only rough 
estimates. As the protein has size of a few nanometers, we take Rq = 10 nm, a = 10 nm and 
Ly = Lz = 30 nm. Temperature is taken to be 300 K . We take the diffusion coefficient of 
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FIG. 1: Plot of the diffusion coefficient of the macromolecule against the longest relaxation time (r). Dp is 
the diffusion coefficient outside the gel. Do and D with a — 2 x IQ^^^J/m (for the gel) and Dbrush with 
0- = 2 X lO^^V/m. 

the protein outside the network to be Dp = ksTlQp = lO^^^m^s"^, which gives the value of 
Qp = 4.14 X J s / m?' . Considering a slab of dimensions Lx,Ly and Lz{= Ly), the energy 

required to deform it by a constant amount a in the Z-direction is where Y is the Young's 

modulus. Using our equation for Edis for the same situation gives kLxa^/2. Equating the two, 
we get k = Y. Y for the hydrogel formed by nucleoporins in vitro has been measured to be 2000 
Pa [13]. With this value for k, to produce a deformation with = 10 nm over a length of 20 nm 
requires only an energy of 1.2kJ/mol, which is ~ ksT/l. So the deformation of the gel to produce 
a hole of the required size requires only thermal energy. So then, why does not the particle just 
deform the gel and get inside? The answer must be that to get inside, it will have to break the 



network. It is known that a colloidal particle of size 5 nm is not able to get into the pore 

This means that there must be at least one hydrophobic contact to be broken to create a hole of 

volume (5nm)^. We will take the energy of a hydrophobic contact as bksT. Therefore, to create a 

hole of the size of a protein, (20nm x lOnm x lOnm) , one will have to break 16 hydrophobic contacts 

which will require an energy of SOksT. This rather large energy requirement can be compensated 

by the formation of roughly the same number of hydrophobic contacts between the protein and 

the network. We estimated a by putting the condition that the deformation energy due to the 

two quadratic terms of E^ig have the same value for creation of a deformation of size 10 nm over 

a length of 20 nm. This gave a for the gel agei = 2 x 10^^'^J/m. As no experimental data on 

the relaxation time of the network is available, we took C = 2 Js/m"^, so that the long wave length 

relaxation time r of the gel would be 1 ms. This seems reasonable as the gel is formed due to 

rather weak hydrophobic interactions among FG units. These data correspond to the following 

values for the dimensionless parameters: a = l,k = l,(p = 0.0207, Dp = 1/Cp = 48.3092. With 

these values, we find: The value of Dq, diffusion coefficient of the protein inside the network to be 

6.51 X W-^^m^/s. This means that putting into the network has reduced the diffusion coefficient 

□ 

by roughly a factor of ~ 15, as expected|10l]. The value of r = 2.4 and hence D = 1.56 x 10 ^'^w? j s 
- noise in the gel enhances the diffusion considerably. Putting the length of the NPC to be 50nm, 
this would give a residence time of 1.6 ms for the particle. This estimates a transport rate of 
~ 600 proteins/second which is roughly in agreement with the experiments [8|]. In Fig. [T]we have 
made a plot of diffusion coefficient against the relaxation time (r) of the gel, keeping all other 
parameters fixed which shows that the faster is the relaxation of the gel, faster is the diffusion of 
the macromolecule (protein) inside it. Further, the enhancement factor, r is ~ 2 for r > 1 ms. 

What would happen if the plug were a brush? In a brush, there are no hydrophobic contacts 
and hence there is no network. Any contact would contribute to the elasticity of the plug and 
therefore, if the plug were a brush, the value of a would be much lower. The model is easily 
analyzed in the limit o" — > 0, to find that Dq ^ 0, r — > 1 so that Z) ^ 0. Thus the diffusion 
coefficient of the particle within the brush would be much lower than in the gel. To demonstrate 
this, we have calculated the diffusion coefficient for the case with a = CTgei/lOO and the results are 
given in the Fig. [H and it is seen that the diffusion coefficient is lowered by a factor of 1/35. 

It is also interesting to compare our results with those of Zilman et al 123\. In their model, the 
greater binding energy of the particle to the FG units would enhance the transport, by retaining 
the particles longer in the pore. The same thing would happen within our model too. The two 
models are similar in that the movement of the particle is co-operative, with the contacts at the 
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front being formed at the same time as the contacts at the back are being broken, leading to easy 
motion of the particle within the NPC. However, over approach is more detailed, and is able to 
take care of the limit where a is small, which as we have argued, is very important. 

In summary, we have proposed a model for the transport of a protein through the NPC. We 
find that a reversible gel would lead to larger diffusion coefficients than the polymer brush and 
hence we suggest that the plug of the NPC is a gel. However, final confirmation of this needs more 
experimental work and simulations. 
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